Energy and angular momentum radiated for non head-on 

binary black hole collisions 
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Abstract 

We investigate the possible total radiated energy pro- 
duced by a binary black hole system containing non- 
vanishing total angular momentum. For the scenarios 
considered we find that the total radiated energy does 
not exceed 1%. Additionally we explore the gravitational 
radiation field and the variation of angular momentum 
in the process. 

After the formation of the final black hole, the model 
uses the Robinson- Trautman (RT) spacetimes as back- 
ground. The evolution of perturbed RT geometries is 
carried out numerically. 

PACS: 04.25.Dm, 04.30.-w, 04.30.Db, 02.60.Cb 



1 Introduction 

The advent of powerful detectors capable of directly 
measuring gravitational radiation for the first time has 
motivated numerous investigations of systems likely to 
produce gravitational waves of enough intensity which 
are expected to be observed with these detectors. Among 
those systems, a prime candidate is that composed by a 
couple of inspiraling black holes which eventually merge 
into one releasing considerable amounts of energy via 
gravitational radiation. Clearly, because of the strong 
gravitational fields involved in such process, its complete 
description requires solving Einstein equations in their 
full generality which can only be done through numerical 
methods jl]]. Several groups are combining efforts to nu- 
merically model such system^], [|||, H|[| and although 
significant progress has been (and is being) made in this 
direction (see for instanccjl]]), an accurate and robust 
implementation is still missing. In the mean time, valu- 
able insight can be gained through approximate models 
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of the system. These models serve both as means to gain 
a better understanding of the problem and also to pro- 
vide information that can serve as tests for the numerical 
simulations, where the lack of known solutions renders 
the problem of accuracy assessment more involved. 

Traditional models are obtained through perturba- 
tive approximations where an expansion with respect 
to some appropriately chosen parameter provides a re- 
duced and manageable system of equations. Naturally, 
these approaches provide accurate answers only when 
the perturbative parameter remains small. For instance, 
the traditional post-Newtonian approximation (see for 
instance H |9|), can be safely used when the black holes 
are far enough and the relative speeds involved are much 
smaller than 1. However, as the holes come closer, the 
gravitational fields and relative speeds involved become 
very large and it is clear that this approximation will 
fail to give sensible answers. Resumation techniques |hJ 
and/or effective one-body expansions [jlT|| have been pro- 
posed to extend the regime of validity of this approach. 
However, even when the ambiguities proper of these op- 
tions can be satisfactorily addressed, they will inevitably 
break in the late stages of the merger, which for the case 
of black holes lies in the maximum sensitivity window of 
earth-based detectors [fr2|, 

Another approximation, known as the close limit ap- 
proximation (CLA) (sec for instance fl4|, |l5| ) assumes 
the black holes have already merged and the pertur- 
bation parameter can be considered to be how non- 
Schwarzschild (or non-Kerr) the hole is. Since no-hair 
theorems imply the final fate of the formed black hole 
should be of the Kerr-type, this approach can safely be 
used to describe such epoch (note however, that in cer- 
tain situations, the CLA has been able to produce valu- 
able results at rather earlier times when the holes had 
just merged; however, it is not clear that this will hold 
in generic cases). 

A different approximate model can be obtained by 
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considering Robinson- Trautman spacetimes [fl6| which 
contain purely outgoing gravitational radiation and de- 
cay to leave a Schwarzschild-like horizon jlT], |l8|. 
These spacetimes are thus natural candidates to study 
systems settling to a single black hole. The idea here 
would be to use them as background spacetimes for a 
perturbation treatment^]. For instance, a perturbation 
approach can be constructed where the perturbation pa- 
rameter is the gravitational radiation produced by the 
system This indeed appears as an enticing sugges- 
tion since the system is not expected to produce total 
radiation outputs in excess of few percent of the initial 
mass. 

In the past, this idea was pursued to produce esti- 
mates of the total gravitational radiation) 21 2^] obtain- 
ing excellent agreement with numerical relativity results 
for the head-on collision, i.e. zero angular momentum 
case) 23 . These studies did not allow for the spacctimc 



to have total angular momentum and also the physical 
situations considered produced small enough radiation 
output that the equations governing the dynamic of the 
perturbations were solved in the asymptotic late-time 
linear regime |2l|, |22| . 

In order to consider more interesting scenarios, pertur- 
bations allowing for angular momentum were introduced 
in pif thus paving the way for studies of more generic 
situations. In the present work, we perform such studies 
by numerically solving the dynamical equations intro- 
duced in [^4| in their full generality. Our studies arc 
aimed towards obtaining a deeper understanding of the 
system and also with the hope of producing information 
that can be used as a check of full numerical relativistic 
simulations. 

This work is organized as follows, in section || we 
briefly review the perturbative treatment of Robinson- 
Trautman (RT) spacetimes that allows for angular mo- 
mentum. Section ^ describes the physical parameters 
and how to make the matching between the two eras 
that appear in our model. The numerical implementa- 
tion to solve the equations that govern the evolution of 
the perturbations is described in Section ^| The numer- 
ical results are presented in Section || where the total 
energy, the gravitational radiation field and the evolu- 
tion of angular momentum are shown. We conclude with 
some final comments in Section ^. 



null tetrad (£ a , m a , m a , n a ) where: 

gab T n b = -g ab m a m b = 1 



(1) 



with all other possible scalar products being zero; then, 
the metric can be expressed by 

g a b = 4"6 + n a 4 - m a m b - m a m b . (2) 

In terms of the coordinate system x\ x 2 , x 3 ) = 
(w, r, (( + C), i(C — 0) > where u is a null coordinate and 
r is an affine parameter along the geodesic integral lines 
of the vector £ a , one can express the null tetrad in terms 
of its components by the relations: 
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with i = 0,2,3 and C = | i 3 ? + * £ 3 ) ! an d where the 
components U and X % are given by the following 
expressions: 
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2 General angular momentum 
perturbations of RT spacetimes 

In generic perturbations of RT spacetimes, which 
include angular momentum, have been presented. 

It is convenient to express the geometry in terms of a 



1 Note that this approach would therefore have a background 
spacetime which radiates as opposed to considering perturbations 
off a stationary spacetime. 
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where 

c.c means complex conjugate, 

K v = ~3v9yV - -^3y V B v V + V 2 



(14) 



and where in these equations we are explicitly denot- 
ing the first order terms by introducing the first or- 
der parameter A dependency, and where 8y is the edth 
operator^ |27|, in the GHP notation |28|], of the 
sphere with metric 



ds 2 = -^d(d<; 



(15) 



where P = V(u, £, C)-Po(C> C)i and Po is the value of P 
for the unit sphere. 

The scalar V is given by 



V = V RT + XV X ; 



(16) 



where Vrt is the RT scalar satisfying the Robinson- 
Trautman equation 



5^2 



-3 Mo Vrt = Vrt 3 3 Vkr- Vg r 3^Vr T 3 Vrt, (17) 

Vx is the linear perturbation scalar and 3 is the edth 
operator of the unit sphere. One can then express Ky 
by 



where 
K Vrt = 

and 
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In the above equations we have 00 = cro(u, £, C), W"o 
WoKnCO, *? = *?(«, CO and 

*3 = *2(«,C,C) 



M + A(Mi(«,c,C) + *>(«,C,C)) ; (21) 



where /1 is related to <7o by 
1 

2~i 



(22) 



In order to study the intrinsic fields at future null in- 
finity (scri), it is convenient to consider the leading order 
behavior of Wo , namely 



4!r 



(23) 



where = C-C) and Wi = 0(l/r 2 ). It is interest- 

ing to note that the component "to = %pr- of the Weyl 
tensor is given in this case[B4l by 



*g d 4 IVi 
r 5 dr 4 ' 



(24) 



where the second term is of order 0(l/r 6 ). 
The remaining equations at scri are: 



*o = + ®v RT *i ~ 3M a , (25) 

VRT 
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The objective is to use the solutions of the previous 
equations to model a black hole that just formed after the 
non head-on collision of a previous binary system. In this 
model the idea is to only make use of the information of 
the individual masses and the total angular momentum. 
With all this in mind we have to chose the appropriate 
gauge and free functions. 

In reference it was discussed the gauge freedom of 
these spacetimes. In order not to introduce extra struc- 
ture, we will chose the free functions and gauges that 
make 

Mi = 0; (28) 
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and 

(To = 0. (29) 

Other choices will force extra structure on the model 
beyond masses and angular momentum; as mentioned 
above. 

With this choice the previous equations become 



®°n = 3 



Vrt 



dv RT K 



(30) 



where A is the amplitude of the excitation, Y 2 fi is an 
spherical harmonic, and (C,C) are the coordinates of 
the sphere where the pole is along the y axis. Then in 
terms of the usual spherical harmonics with coordinates 
(C, C)> one nas 




-y 2 ,_ 2 (C,C) + ^2,o(C,C) 



5*2,2 (C, C) I 

(40) 



*o = 3 Vra^o 



and 



( 31 ) 3 Fixing the parameters and the 
Newtonian matching 



-6M Vx = V rt B Vrt Q Vrt K Vx - 18 MqVrt 



V\ 
Vrt 



(32) 

It is interesting to note that since the RT equation can 
also be expressed as 

-6M 0T ^ = ^Vrt^Vrt k v rt , (33) 
Vrt 

equation (|3^) can be written in the following way 

- 6M V\ = Vrt^Vrt^Vrt k v x 

+ 3 VxEv RT 9y„ T K Vrt . (34) 

Given initial conditions for the functions Vrt, ^q, ^1 
and V x equations @, ©, @ and © provide for 
the respective evolutions. 

The radiation content of the spacetime is easily de- 
scribed by the time derivative of the Bondi shear. Let 
a B denote the Bondi shear and ug the Bondi time; then 
for any section of scri u = constant one can express the 
radiation content from the relation |E9| 



Y 



8ub V 
therefore to first order one has 



A 



d 2 Vx VxB 2 Vrt\ 



V 2 

RT 



-) 



(35) 



(36) 



du B Vrt ' ' V Vrt 
The radiation flux at the retarded time u is given by 
1 



F B (u) 



0a% B-al ds2 _ 



47T Jg u 8ub dus 



(37) 



Let us recall that in this case the timelike component 
of the Bondi momentum at the section u = constant 
can be calculated from 



K = ±[^ds 2 . 



4tt / V 3 



(38) 



At the retarded time it, we say that V represents a 
quadrupole excitation along the y axis if 



V = l + AY 2t0 (C,t')\ 



(39) 



The physical system can be characterized by two stages. 
During the first stage two black holes are moving to- 
wards each other with some orbital angular momentum. 
At some moment they collide and form a single black 
hole; which settles down in the asymptotic future to a 
stationary Kerr geometry. 

During the first stage of evolution we describe the 
gravitational radiation with the quadrupole formula, 
where the dynamics is worked out from the Newtonian 
framework. The whole motion is contained in a plane; 
therefore two variables are sufficient to describe the or- 
bits. The two integrals of motion, namely energy E and 
angular momentum j, allow to solve the Newtonian sys- 
tem. 

The initial data is assumed to be given in the (x, y) 
plane; with zero total momentum and such that the or- 
bital angular momentum is along the positive z direction 

The initial velocities can be thought to have compo- 
nents along the x axis only. Let us call Rq the initial 
impact parameter, at some initial relative distance tq. 
In this way, using the relative velocity Vq, one has 



E 



j = fiR v ] 
is the reduced mass. 



(41) 



(42) 



where a = —, 

For this kind of motion the quadrupole formula pre- 
dicts that the power of gravitational radiation is given 
by 



F Q {r) = — {Gm 1 m 2 y 



fir"* 



E 



G 



r 



fi 2 r 6 



. (43) 

which generalizes the analog equation appearing in ref. 
II- 

In order to continue the dynamical description after 
the collision of the two black holes, we need to have a 
merging condition in the Newtonian framework. In ref- 
erences |2l| and ||^] we have succeeded in estimating the 
total energy radiated in the head on black hole collision 
using the following criteria. When a black hole, which 
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mass at infinity is win, is brought to a distance m of 
another black hole of asymptotic mass m,2, its physical 
mass, for the stationary situation, is changed ]3C|| to 



ml = run 



2r 12 : 



(44) 



and similarly for the other mass. 

Applying the arguments of rcf. one concludes that 
the merging condition should be taken when the separa- 
tion distance has the value 



r c = 2 (ml + ml) . 



(45) 



It is convenient to introduce the relative mass param- 
eter a and the reference mass m, so that mn = m and 
771^2 = am. Then the merging condition p2| gives 



(mil + m l2 ) 




(46) 



Let us now consider the initial data used in rcf. Q for 
the grazing collisions of black holes; namely: x\ = 5m, 
x 2 = —5m, yi = m,y 2 — — m, v xl — —0.5 and v x2 = 0.5: 
with m = 1 and a = 1. This data corresponds to a 
hyperbolic Newtonian trajectory; with initial separation 
distance r = 10.198, and critical merging radius r c = 
4.449. 

Our strategy is to follow closely the model used in 
[ pTj , |22| ; but we also want to compare our work with Q . 
However, the initial data of Q involves two black holes 
with half the speed of light each; therefore in adapting 
this initial data to our model we must take into account 
relativistic effects. In references [^lj and there was 
no need for these concerns because the initial data was 
not relativistic. 

Since the initial velocities are relativistic, the relative 
initial velocity vq is calculated from the expression 



Vx2 — Vxl 

1 - v x iv x2 



0.8; 



(47) 



where we are using geometric units in which the veloc- 
ity of light and the gravitational constant have the unit 
value. 

The energy radiated during the falling phase, calcu- 
lated from @ is E N = 0.00346(mi + m 2 ). 

In order to match the Newtonian stage with the black 
hole RT perturbed model we need also to set the total 
initial mass and initial angular momentum for the RT 
stage. 

In previous work we have matched the Newtonian 
mass M]y ew = mn + m.i 2 to the initial mass M; this 
already takes into account the field relativistic first or- 
der correction in terms of its physical masses ml and m2, 
since, recalling equation (21) of rcf. (22), the initial mass 
would be ml + m2 — m \]™ i2 ■ The system discussed in 
ref. JIH had zero initial velocity, and therefore there was 
no need to take into account any other effect. Instead in 



our case we should take into account speed relativistic 
corrections. 

Then, since the initial data is relativistic, the initial 
mass and angular momentum are calculated from 



M 



and 



-tJi 



m,i\v x \ 



+ 



m i2 



= 2.309 



m t2 v x2 
U2 , = 1.155. 



(48) 



(49) 



-:.:2 



It should be emphasized that the relation between the 
initial RT stage mass and angular momentum M and J 
with the Newtonian mass and angular momentum M^ ew 



and j is just the Lorentzian factor 7 = 



1.1547. 



In the last section we will comment on the incidence of 
this factor on our results. 

Let us observe that the relation between the relativis- 
tic angular momentum and total mass is jj? = 0.217. 

Since the angular momentum is small, it can be 
treated as a perturbation. At the moment of the col- 
lapse, we can consider a quadrupole excitation along the 
y axis with amplitude A, as described above. Let us 
take A = A m + Aj ; where Aj is the contribution com- 
ing from the appearance of the angular momentum j. 
The matching condition is given by the equation (See 
equations (||), @ and ©.) 



F B =F Q (r c ); 



from which one obtains in this case 



A rn = 0.057 



and 



Aj = 0.028. 



(50) 



(51) 



(52) 



Therefore, to be explicit, the initial data for Vrt and V\ 
are 



Vrt = l~-f- 

x (y|^_ 2 (CO + ^o(C, + ^1^2,2 (CO) ; (53) 



and 



A Vx = — 



A ± 
2 




i((,0 + Y 2 ,o(CX) 



2^, 2 (C,C) 



(54) 



The constant M is determined from ( |3q) , and the 
condition that initially the mass is given by ([48|); which 
sets M = 2.302. 
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The orbital angular momentum is taken into account 
in the initial data for W". It is convenient to express this 
initial data in terms of the auxiliary field g, given by 



and corrector 



V 2 



Vrt 



(55) 



Let us note that then 



so that in the stationary case one has 9 2 <? = and g = 0. 
We take g = g and 



9 = .9o^i,o(CiC); 



(56) 



where go is related to the angular momentum in the z 
direction by 



3o 



-4.094. 



(57) 



In reference |3l| it was also considered a similar case of 
a binary system with orbital angular momentum. Their 
initial data was: x\ = 0, X2 = 0, y\ = to, yi = —to, 
V x i = —0.8 and v X 2 = 0.894; with m =1.5 and a = |. 
From the Newtonian point of view this data corresponds 
to an elliptic motion; but with maximum and minimum 
radius that arc smaller than the corresponding critical 
merging radius. For this reason, we can not compare 
this case with our model. 



4 Numerical Implementation 

Accurate numerical evolution of a fourth order parabolic 
equation, such as (|l7|), by means of an explicit finite 
difference scheme is a challenge because the CFL|32| 
condition requires that the time step Am scale as the 
fourth power of the spatial grid size. Nevertheless, we 
constructed a set of algorithms to solve these equations 
using second order accurate finite difference approxima- 
tions (following |53|). The numerical treatment of the 
eth operator has been thoroughly described in Q . This 
work presented a clean way to deal with derivative op- 
erators on the sphere by covering it with two coordinate 
patches and dealing with spin weighted quantities. Thus, 
it is ideally suited for our present purposes. The numer- 
ical grid on each patch is defined by £y = qi + ipj where 
q hPi = -l-2A A + {i-l)A A (with A A = 2/(N A -5)). 
The angular derivatives are discretized by centered sec- 
ond order finite difference approximations and informa- 
tion between patches is obtained through fourth order 
accurate interpolations. (For a detailed description of 
this approach sceQ). 

The integration in time is based upon a three time 
level Adams-Bashford [|D scheme with predictor {£) 
given by 



P{u + Au) = T{u) + —d u [6T{u) - T{u - Au)}, (58) 



T(u + Au) = T(u) + —dul^iu) + T{u + Au)} 

+ 0(Au 3 ). (59) 

Where T stands for Vrt or V\ and the d u terms 
are to right hand sides of equations d27]). Addition- 
ally, we implemented the iterative Cranck-Nicholson 
algorithm |3^] and observed that the results obtained 
with both implementations agree. Since the evolution 
equation for \1/q is linear, its numerical integration is 
straightforwardly done by centered second order differ- 
ences at the level (u + Au/2). 

The second order convergence of numerical solutions 
was confirmed in the perturbative regime using solu- 
tions of the linearized equation and second order self- 
convergence of the solutions was confirmed in the non- 
linear regime. 



5 Results 

Variations of total energy 

The total Bondi energy-momentum vector at any RT 
retarded time u can be calculated from the expression 



P * = -T- I nCX) (*°B2 + <rB*B) dS 2 

J s 



(60) 



where S is the section determined by u = constant, 

c + c c - c cc - 1 



i + cc'*(i + ccT i + cc 



(61) 



and the subscript B is used to emphasize that the quan- 
tities are evaluated with respect to a Bondi frame. The 
mass M at this section S is then given by 



(62) 



where the indices are raised and lowered by the 
Lorentzian flat metric r/ a b at scri p6| . 

Let us note that the relations between the Bondi quan- 
tities and the intrinsic ones are 



and 



a 



(63) 



(64) 



therefore in our gauge one has as — 0, at each RT sec- 
tion. 

The gravitational energy radiation flux is calculated 
from the Bondi time derivative of the supcrmomcntum 
frpTll; namely 



das das 



duB 8ub dus 



(65) 
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If one instead considers the time change with respect to 
the RT time, it is convenient to have in mind that for 
any function / one has 



df_ 

dun 



l_df_ 

V du 



(66) 



The so called news function can be expressed in 

OU B 1 

terms of the perturbed RT fields by 



d<J B 
du B 




(67) 



To calculate the total energy radiated, one could then 
numerically evaluate the gravitational energy radiation 
flux of equation (^) at different times and sum along 
all the elapsed time. However, it is more accurate to 
numerically evaluate the initial mass and subtract the 
final mass. This is due to the fact that the RT spacetime 
is known to converge asymptotically to the Schwarzschild 
one; more specifically, one knows that linin^oo Vrt = 1; 
and similarly one can see that lim u _>oo V\ = 0. 

Using this procedure, and a resolution of n = 32 
points for half a meridian of the sphere (approximately 
TV = 1600 points for the whole sphere), the energy ra- 
diated Ert during the RT stage is found to be Ert = 
0.0034M . 

Then, since in units of Mo, the energy radiated in the 
first stage is Ejy = 0.0030M , the total energy radiated 
in the whole process is E = 0.0064M . 

Gravitational radiation field 

The numerical calculation of the evolution of the grav- 
itational radiation field ^4, is shown in figure [j]. 

It can be seen that although the total energy radiated 
is rather small, the amplitude of the gravitational ra- 
diation field can be large. In other words, this model 
describes a noticeable burst. 

This is interesting since ^4 is precisely what gravita- 
tional wave detectors will measure. 

Variations of angular momentum 

When dealing with the notion of angular momentum 
one is faced with the fact that there are several inequiv- 
alent definitions of angular momentum; which are not 
tightly related with the notion of intrinsic angular mo- 
mentum, with the exception of [fj8f . An appropriate def- 
inition of intrinsic angular momentum involves the selec- 
tion of unique sections J38| of future null infinity where 
the quantity is to be calculated. Since in our case we 
are taking the angular momentum as a perturbation pa- 
rameter of the RT geometries, it is not essential to con- 
sider these refinements in our model. And also, since the 
RT spacetimes provide with a geometric unique family 
of sections of future null infinity, namely the sections 
u = constant, it is natural to use them to calculate the 
angular momentum. 

Then, instead of describing the variation of the intrin- 
sic angular momentum we describe the variation of the 



Re(T 4 y) 



-0.1 



-0.2 

0.4 
0.3 
0.2 
0.1 

-0.1 
-0.2 







0° 

30° 

60° 
90° 


/ 










/ -0.1 






%/ -0.2 






0.001 



0.2 0.4 0.6 0.8 



1.2 1.4 1.6 1.i 



Re(W 4 x) 




U 
30° 
60° 



0.001 



0.2 0.4 0.6 0.! 



1 1.2 1.4 1.6 1.i 



Figure 1: Evolution of ^4 as a function of time. The 
subscript x and y refer to the x and y axis of the 
stereographic coordinate £ = x + iy of the north- 
ern hemisphere. Each of the four curves refer to 
detectors at 0°, 30°, 60° and 90°; measured from 
the north pole of the sphere. The insets show the 
detail of the rapid time variations at the beginning. 



RT-angular momentum vector given by 



[ ±M k *° m dS 



(68) 



where Srt are the sections determined by u = 
constant, k = 1,2,3, so that t k are the spacelike com- 
ponents of £ a ; and where the Bondi component ^S B1 , of 
the Weyl tensor, is related to the RT Weyl component 
*S by 

*§ 
W 

Figure ^| shows a very small and smooth variation of 
the angular momentum; which is more related to the 
time variation of the RT geometry than to the radiation 
of angular momentum; as can be seen from the nature 
of equation (^Tj). 



(69) 
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Figure 2: Evolution of the relative variation of an- 
gular momentum J as a function of time. 

6 Final comments 

The estimate of the total gravitational energy radiated 
in the non head-on collision of two black holes is less than 
one percent according to our model. This is on the low 
end value of existing numerical results (which are still 
being refined) and comparable with those obtained with 
the close limit approximation |59j. It is however worth 
mentioning that our calculations predict an amplitude 
of the gravitational waves that could be important from 
the observational point of view. Additionally, it predicts 
a rather narrow time duration of the 'burst'. Since the 
frequency at which this happens; its amplitude and the 
time duration can be used as preliminary information for 
constructing data analysis 'filters'^, |4l] [, knowledge of 
these is of importance while templates from full numer- 
ical simulations are not available. We will carry out a 
detailed study of these in the future. 

At first sight there is too large a difference between 
the numerical results in j^] and the estimates obtained 
here. However, one should keep in mind several points 
which (by themselves) can account for this difference. 
Comments are in place for both the results in ^ and 
those presented here: 

First, the results in ^ were obtained under a rather 
coarse resolution and the final numbers need to be re- 
fined (there are intense efforts worldwide in this direc- 
tion). Additionally, the obtained value of radiated en- 
ergy come from comparisons of initial ADM mass and 
apparent horizon masses. Masses obtained from appar- 
ent horizon calculations, in dynamical regimes, are only 
an approximation of the mass the black hole and hence 
results obtained through this method can have a signifi- 
cant systematic error. Furthermore it is important to re- 
member that the initial data of [|| is given at a spacelike 
hypersurface of the spacetime; while, by the nature of the 
RT spacetime, our initial data is given on a characteris- 
tic surface. This implies an essential difference since in 



the spacclikc hypersurface initial data one has incoming 
and outgoing radiation contributing to the total energy 
radiated; while in the RT characteristic problem that we 
solve, we do not consider any incoming initial radiation 
in our calculations. So, for the relativistic regime, one 
expects the time slice initial data to radiate more energy 
than the characteristic case. 

As far as the present model is concerned: in previous 
works we have obtained very good agreements between 
our estimates and the mature numerical calculations of 
exact geometries; as one can check in the family of sys- 
tems depicted in figures 3 and 4 of reference [Q. The 
fact that the model represented so well such a variety 
of initial data and different mass ratios, motivated us to 
apply the same techniques to the case in which angular 
momentum was involved. However, since there are just 
two (still being refined) numerical calculations available 
to compare with, it is difficult to obtain much informa- 
tion at this stage. 

As we have mentioned before our aim is to follow as 
close as possible the model used in [^l] and ; however 
the fact that the initial conditions considered are rela- 
tivistic posse several questions. For example, in match- 
ing the masses of the two stages we have used equation 
((f8[); so one might ask how would the results change if 
the relativistic 7 factor was not used? This would cer- 
tainly decrease the initial M and Mq by about 13%; but 
En and Ert would not change, so that the ration be- 
tween the total energy radiated versus the changed Mq 
would increase by the factor 7, that is about 15%. 

Similarly, had we not used the relativistic correction 
on the angular momentum J, its initial value would have 
decreased by about 13%; but the relative variation shown 
in figure |^ would not have changed. 

In any case, we think that by applying our model to 
this relativistic initial data, we are pushing the model 
to the boundary of its validity. Comparisons with future 
numerical simulations will shed light on this and indicate 
how far this model can be pushed. 

The radiation of angular momentum seems to be neg- 
ligible with these initial conditions. In order to consider 
higher values of the initial angular momentum, we would 
need to deal with other background geometries, as for 
example twisting algebraically special spacetimes. Re- 
garding the smooth monotonic variation of it one can 
infer that, for this small initial angular momentum data, 
its behavior is driven by the exponential asymptotic be- 
havior of the RT background geometry. There are not 
complicated initial variations in '5 J, that for example do 
appear in ^4, as seen in figure [j]. 

When describing a concrete physical situation with 
these spacetimes, one is supposed to choose the gauge 
and fix the free functions in order to make the best rep- 
resentation of the system. It is somehow striking that 
the choice of the frame in first order has physical sig- 
nificance, and it is not pure gauge as one is accustomed 
in the studies of linearized gravity around Minkowski 
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spacctime. 
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